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Abstract 



We present new solutions in terms of elementary functions of the multi-component non- 
linear Schrodinger equations and known solutions of the Davey-Stewartson equations such as 
multi-soliton, breather, dromion and lump solutions. These solutions are given in a simple 
dctcrminantal form and are obtained as limiting cases in suitable degenerations of previously 
derived algebro-geometric solutions. In particular we present for the first time breather and 
rational breather solutions of the multi-component nonlinear Schrodinger equations. 

1 Introduction 

One of the significant advances in mathematical physics at the end of the 19th century has been 
the discovery by Gardner, Greene, Kruskal and Miura |18] of the applicability of the Inverse 
Scattering Transform (1ST) to the Korteweg-de Vries equation, and the construction of multi- 
soliton solutions. The most important physical property of solitons is that they are localized 
wave packets which survive collisions with other solitons without change of shape. For a guide 
to the vast literature on solitons, see for instance |3H HOj . Existence of soliton solutions to the 
nonlinear Schrodinger equation (NLS) 



where p = ±1, was proved by Zakharov and Shabat [12] using a modification of the 1ST. The NLS 
equation is a famous nonlinear dispersive partial differential equation with many applications, 
e.g. in hydrodynamics (deep water waves), plasma physics and nonlinear fiber optics. The N- 
soliton solutions to both the self-focusing NLS equation (p = 1), as well as the defocusing NLS 
equation {p = —1), can also be computed by Darboux transformations [28J, Hirota's bilinear 
method (see e.g. [20, 34, 9j) or Wronskian techniques (see |17 [ I3"0" l fTB]). Hirota's method relies 
on a transformation of the underlying equation to a bilinear equation. The resulting multi-soliton 
solutions are expressed in the form of polynomials in exponential functions. Wronskian techniques 
formulate the iV-soliton solutions in terms of the Wronskian determinant of N functions. This 
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method allows a straightforward direct check that the obtained solutions satisfy the equation 
since differentiation of a Wronskian is simple. On the other hand, multi-soliton solutions of 
(1.1) can be directly derived from algebro-geometric solutions when the associated hyperelliptic 
Riemann surface degenerates into a Riemann surface of genus zero, see for instance [7]. 

In the present paper, we construct solutions in terms of elementary functions of two gener- 
alizations of the NLS equation (1.1): the multi-component NLS equation (n-NLS), where the 
number of dependent variables is increased, and the Davey-Stewartson equation (DS), an in- 
tegrable generalization to 2 + 1 dimensions. The solutions of n-NLS and DS presented in this 
paper are obtained by degenerating algebro-geometric solutions, previously investigated by the 
author in [22j using Fay's identity [29J. This method for finding solutions in terms of elementary 
functions has not been applied to n-NLS and DS so far. It provides a unified approach to various 
solutions of n-NLS and DS expressed in terms of a simple determinantal form, and allows to 
present new solutions to the multi-component NLS equation in terms of elementary functions. 

One way to generalize the NLS equation is to increase the number of dependent variables in 
(1.1). This leads to the multi-component nonlinear Schrodinger equation 

1 it + + 2 (E s ^I 2 ) *s = °> i = L • • ■ > n ' (i- 2 ) 

denoted by n-NLS s , where s = (s\, . . . ,s n ), = ±1. Here vpj(x,t) are complex valued func- 
tions of the real variables x and t. The case n = 1 corresponds to the NLS equation. The 
two-component NLS equation (n = 2) is relevant in the study of electromagnetic waves in 
optical media in which the electric field has two nontrivial components. Integrability of the 
two-component NLS equation in the case s = (1,1) was first established by Manakov |25j . In 
optical fibers, for arbitrary n > 2, the components tpj in (1.2) correspond to components of the 
electric field transverse to the direction of wave propagation. These components of the transverse 
field form a basis of the polarization states. Integrability for the multi-component case with any 
n > 2 and Sk = ±1 was established in [38] . Multi-soliton solutions of (1.2) were considered in a 
series of papers, see for instance [25j (35j ESI [23j [T|. 

In this paper, we present a family of dark and bright multi-solitons, breather and rational 
breather solutions to the multi-component NLS equation. This appears to be the first time that 
breathers and rational breathers are given for the multi-component case. The notion of a dark 
soliton refers to the fact that the solution tends asymptotically to a non-zero constant, i.e., it 
describes a darkening on a bright background, whereas the bright soliton is a localized bright spot 
being described by a solution that tends asymptotically to zero. The name 'breather' reflects 
the behavior of the profile which is periodic in time or space and localized in space or time. 
It is remarkable that degenerations of algebro-geometric solutions to the multi-component NLS 
equation lead to the breather solutions, well known in the context of the one-component case as 
the soliton on a finite background [3J (breather periodic in space) , the Ma breather [24J (breather 
periodic in time) and the rational breather |33j . In the NLS framework, these solutions have 
been suggested as models for a class of extreme, freak or rogue wave events (see e.g. |19| l32l S]). 
A family of rational solutions to the focusing NLS equation was constructed in |13| and was 
rediscovered recently in [IT] via Wronskian techniques. Here we give for the first time a family 
of breather and rational breather solutions of the multi-component NLS equation. For the one 
component case, our solutions consist of the well known breather and Peregrine breather of the 
focusing NLS equation. For the multi-component case, we find new profiles of breathers and 
rational breathers which do not exist in the scalar case. 
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Another way to generalize the NLS equation is to increase the number of spatial dimensions 
to two. This leads to the DS equations, 



i ik + i>xx - a 2 tp yy + 2 ($ + p I V| 2 ) V = 0, 



yy 

$xx + a 2 $ ra + 2/>|V|L = 0, (1.3) 



where a = i, 1 and p = ±1; if)(x,y,t) and $(x,y,t) are functions of the real variables x, y and 
i, the latter being real valued and the former being complex valued. In what follows, DS1 P 



corresponds to the case a = i, and DS2 P to a = 1. The DS equation (1.3) was introduced 
in |12| to describe the evolution of a three-dimensional wave package on water of finite depth. 
Complete integrability of the equation was shown in (5]. A main feature of equations in 1 + 1 
dimensions is the existence of soliton solutions which are localized in one dimension. Solutions 
of the 2 + 1 dimensional integrable equations which are localized only in one dimension (plane 
solitons) were constructed in [21 [6]. Moreover, various recurrent solutions (the growing-and- 
decaying mode, breather and rational growing-and-decaying mode solutions) were investigated 
in |40j . The spectral theory of soliton type solutions to the DS1 equation (called dromions) with 
exponential fall off in all directions on the plane, and their connection with the initial-boundary 
value problem, have been studied by different methods in a series of papers |8j [151 H22 EZ]- The 
lump solution (a rational non-singular solution) to the DS2 - equation was discovered in p]. 

Here we present a family of dark multi-soliton solutions to the DS1 and DS2 + equations, as 
well as a family of bright multi-solitons for the DS1 and DS2 - equations, obtained by degener- 
ating algebro-geometric solutions. Moreover, a class of breather and rational breather solutions 
of the DS1 equation is given. These solutions have a very similar appearance to those in 1 + 1 
dimensions. In this paper it is shown how the simplest solutions, the dromion and the lump 
solutions can be derived from algebro-geometric solutions. 

The paper is organized as follows: Section 2 contains various facts from the theory of theta 
functions and identities due to Fay. These identities were used to construct algebro-geometric 
solutions of n-NLS and DS equations in [22J, and will be needed for the degeneration of the 
underlying theta-functional solutions. Section 3 provides technical tools dealing with the degen- 
eration of Riemann surfaces. We present a method which allows to degenerate algebro-geometric 
solutions associated to an arbitrary Riemann surface that can be applied to general integrable 
equations. In Section 4 solutions in terms of elementary functions to the complexified n-NLS 
equation are derived by degenerating algebro-geometric solutions; for an appropriate choice of 
the parameters one gets multi-solitonic solutions, and for the first time breather and rational 



breather solutions to the multi-component NLS equation (1.2). In Section 5 a similar program is 
carried out for the DS equations; well known solutions such as multi-solitons, dromion or lump 
are rediscovered from an algebro-geometric approach. 

2 Theta functions and Fay's identity 



Solutions of equations (1.2) and ( |1.3[ ) in terms of the multi-dimensional theta function were dis- 
cussed in [22J. In this section we recall some facts from the construction of these solutions which 
will be used in the following to get particular solutions as limiting cases of algebro-geometric 
solutions. 
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2.1 Theta functions 



Let lZ g be a compact Riemann surface of genus g > 0. Denote by {-Aj,Bj}^ =1 a canonical 
homology basis, and by {wj}| =1 the dual basis of holomorphic differentials normalized via 

/ u>j = 2m6 kt j k,j = l,...,g. (2.1) 

The matrix B of i3-periods of the normalized holomorphic differentials with entries (B)fcj = 
Jg fe ujj is symmetric and has a negative definite real part. The theta function with (half integer) 
characteristic 5 = [5', 6"] is defined by 

6[<J](z|B) = exp{i(B(m + (5 , ),m + (5 / ) + (m + y,z + 2i7ry / )} (2.2) 

for any z G C 9 ; here 5', 5" G {0, are the vectors of characteristic and (., .) denotes the 
scalar product (u, v) = Yli u i v i f° r an Y u > v £ C 9 . The theta function 0[<5](z) is even if the 
characteristic 5 is even i.e, 4(<5',<5") is even, and odd if the characteristic 5 is odd, i.e., 4 {5', 8") 
is odd. An even characteristic is called non-singular if @[<f](0) ^ 0, and an odd characteristic is 
called non-singular if the gradient V@[£](0) is non-zero. 



2.2 Corollaries of Fay's identity 

Let us first introduce some notation. Let k a denote a local parameter near a G lZ g . Consider the 
following expansion of the normalized holomorphic differentials Uj near o, 

<"j(p) = (Ya,j + W a ,j k a {p) + o (k a (p))) dk a {p), (2.3) 

where p lies in a neighbourhood of a, and V a j, W a ,j 6 C Let us denote by D a the operator of 
directional derivative along the vector V a = (V a ,i, • • • , Va^Y'- 

9 

D a F(z) = J2d Zj F(z) V aJ = (VF(z), V«) , (2.4) 
i=i 

where i 7 : C 5 — > C is an arbitrary function, and denote by D' a the operator of directional 
derivative along the vector W a = (Wa,l 5 • • • > Wa^) 4 . 

Now let <5 be a non-singular odd characteristic. For any z G C 9 and any distinct points 
a, & G 7£ g , the following two versions of Fay's identity |14j hold (see |29| and |22| ) 

e(z + r)9(z-r) 
6(zp 



ZWn6(z) = gi + ^^T" ' ( 2 - 5 ) 



^ ln + ^ 111 + ( Da ln " Kl ) 2 + 2 D ' ln 0(z) + K > = °' (2 - 6) 

where the scalars qi, Ki for i = 1, 2 depend on the points a, 6 and are given by 

g 1 (a,6) = J D a L» 6 lnG[<5](r), (2.7) 

i? a e[fl(o)A,e[<5](o) 

ft(a,6) = , (2.8) 
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K 2 (a,b) = -D' a \n@(r)-Dl\n(@(r)@(0)) - [o a In 9(r) - K x {a, b)^ . (2.10) 

Here we used the notation r = J u where oj = (u)\, . . . , ojg) 1, is the vector of the normalized 
holomorphic differentials. 



2.3 Integral representation of q 2 (a,b) and K^a^) 



Quantities q2{a, b) and K\(a, b) defined in (2.8) and (2.9) respectively, admit integral representa- 
tion which will be more convenient for our purposes. These integral representations follow from 
the fact that meromorphic differentials normalized by the condition of vanishing ^.-periods can 
be expressed in terms of theta functions. 

Let a,b G lZ g be two distinct points connected by a contour which does not intersect A and 
i3-cycles. Hence we can define the normalized meromorphic differential of the third kind £lb-a 
which has residue 1 at b and residue —1 at a. Now let a G lZ g , and N G N with N > 1. The 
normalized meromorphic differential of the second kind has only one singularity at the point 
a and is of the form 



1 



+ 0(1) )dka(p), p£K g , 



(2.11) 



Mp) n 

where k a is a local parameter in a neighbourhood of a. 

Proposition 2.1. Let a,b G lZ g be distinct points. Denote by k a and kb local parameters in a 
neighbourhood of a and b respectively. The quantities q2(a,b) and K\{a,b) defined in (2.8) and 
(2.9) respectively admit the following integral representations: 

q2(a,b) -- 



lim 

a— >a 



k a (a) k b (b) 



-i 



exp 




b— a 



(P) 



(2.12) 



where the integration contour does not cross any cycle of canonical basis, and 



K x {a,b) 



lim 

a— >a 



V a 2 \p) + 



k a (a) 



(2.13) 



where c is an arbitrary point on 1Z„. 



Proof of (2.12) can be found in |22| . where similar statements lead to (2.13). 



3 Uniformization map and degenerate Riemann surfaces 

It is well known that solutions in terms of theta functions are almost periodic due to the pe- 
riodicity properties of the theta functions. In the limit when the Riemann surface degenerates 
to a surface of genus zero, periods of the surface diverge, and the theta series breaks down to 
elementary functions. Whereas this procedure is well-known in the case of a hyperelliptic surface, 
i.e., a two-sheeted branched covering of the Riemann sphere, where such a degeneration consists 
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in colliding branch points pairwise, it has not been applied so far to thet a- functional solutions 
on non-hyperelliptic surfaces. 

We present here a method to treat this case based on the uniformization theorem for Riemann 
surfaces. In particular, we show that the theta function tends to a finite sum of exponentials in 
the limit when the arithmetic genus of the associated Riemann surface drops to zero, and give 
explicitly the constants (2.7)-(2.10) in this limit. As illustrated in Section 4 and 5, particular 
solutions of n-NLS and DS such as multi-solitons, well known in the theory of soliton equations, 
arise from such degenerations of algebro-geometric solutions. 



3.1 Degeneration to genus zero 

Let us first recall some techniques used for degenerating Riemann surfaces (see [14J for more 
details). There exist basically two ways for degenerating a Riemann surface by pinching a cycle: 
a cycle homologous to zero in the first case, and a cycle non-homologous to zero in the second 
case. The first degeneration leads to two Riemann surfaces whose genera add up to the genus of 
the pinched surface, whereas the limiting situation for the second degeneration is one Riemann 
surface of genus g—1 with two points identified, g being the genus of the non-degenerated surface. 
In both cases, locally one can identify the pinched region to a hyperboloid 

y 2 = x 2 -e, (3.1) 

where e > is a small parameter, such that the vanishing cycle coincides with the homology 
class of a closed contour around the cut [—y/e, \/e\ in the x-plane. In what follows, we deal with 
the degeneration of the second type and make consecutive pinches until the surface degenerates 
to genus zero. 

To degenerate the Riemann surface 7l g of genus g into a Riemann surface TZq of genus zero, 
we pinch all ^.j-cycles into double points. After desingularization one gets TZo, and each double 
point corresponds to two different points on TZq, denoted by it, and Vi for i = 1, . .. ,g. In this 
limit, holomorphic normalized differentials Ui become normalized differentials of the third kind 
with poles at Ui and vi. Note that the normalized differential of the second kind 0>a with a 
pole of order N > 1 at a remains a differential of the second kind with the same order of the pole 
after degeneration to genus zero. We keep the same notation for the differential of the second 
kind on the degenerated surface. 

The compact Riemann surface TZq of genus zero is conformally equivalent to the Riemann 
sphere with the coordinate w. This mapping between TZo an d the w-sphere is called the uni- 
formization map and we denote it by w(p) = w for any p £ TZq. Therefore, in what follows we 
let TZq stand also for the Riemann sphere with the coordinate w. 

Meromorphic differentials on TZq can be constructed using the fact that in genus zero, such 
differentials are entirely defined by their behaviors near their singularities. This leads to the 
following third and second kind differentials on TZq: 



Differentials of the third kind: 



n Vi - Ui = ( 1 1 ) dw. (3.2) 

w - w Vi w - w Ui 



• Differentials of the second kind: 
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1 

k' a (w a 



(3.3) 



(W - W a ) 2 ' 

where k a is a local parameter in a neighbourhood of w a £ TZq and the prime denotes the d erivative 



(2) 

with respect to the argument. This is the differential on TZq, obtained from Qa (2.11) defined 
on TZg, in the limit as the surface 7Z g degenerates to TZq. The factor (^(w a )) _1 ensures that the 

(2) 

biresidue of Qa with respect to the local parameter k a is 1 as before the degeneration. 



3.2 Degenerate theta function 

To study the theta function with zero characteristic in the limit when the genus tends to zero, 
let us first analyse the behavior of the matrix B of ,8-periods of the normalized holomorphic 
differentials. Since holomorphic normalized differentials uii become differentials of the third kind 
with poles at it, and «j, for a small parameter e > 0, elements (B)jfc of the matrix B have the 
following behavior 



) ik = / &v k -u k + 0(e), i / k, 

J Ui 

) fcfc = lne + 0(l). 



(3.4) 



Therefore, the real parts of diagonal terms of the Riemann matrix tend to — oo when e tends to 
zero, that is when the Riemann surface degenerates into the Riemann surface TZq. It follows that 



the theta function (2.2) with zero characteristic tends to one, since only the term corresponding 
to the vector m = in the series may give a non-zero contribution. 



To get non constant solutions of ( 1.2 ) and ( 1.3 ) after the degeneration of the Riemann surface, 
let us write the argument of the theta-function in the form Z — D, where D is a vector with 
components D k = (1/2) (B)^fc + d^, for some d k £ C independent of e. Hence for any Z £ C 9 
one gets 



lim6(Z 



me{0,l} 9 \l<i<k<g 
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E 

fc=i 



ik mim k +) m k (Z k - dk) 



(3.5) 



Here we use the same notation for the quantities 



lik 



on the degenerated surface. The expression 



in the right hand side of (3.5) can be put into a determinantal form (see Proposition 3.1 ) which 
will be used in the whole paper. This determinantal form can be obtained from the following 



representation of the components (B)jfc after degeneration, obtained from (3.2) and (3.4) 



Hence, following [27] one gets 
Proposition 3.1. For any z 6 



)ik 



In 



w. 



w 



W„ 



IV 



■t>k 



the following holds 



(3.6) 



E exp 

me{0,l} 9 
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(b)^ ///,////, + 2^ '"/ ?k 



l<i<k<g 



k=l 



det(T), 



(3.7) 
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where T is a g x g matrix with entries 

(T) ik = S i;k + Wvi ~ Wui e *("+*>. (3.8) 

w„ 4 - 



3.3 Degenerate constants 



The next step is to give explicitly the quantities (independent of the vector z) appearing in (2.5 ) 



and (2.6), i.e., V a , W a ,r, q 2 , etc, after the degeneration to genus zero. We use the same notation 



for these quantities on the degenerated surface. For any distinct points a, b E IZq, it follows from 



(|2_3j) and (|3_2j) that 

Va,k- 

W ajk : 



k' a (Wa) \W a -W Vk W a -W Uk 



+ 



K(Wa) 2 \ (w a -W Vk ) 2 {w a -W Uk ) 



Kb" a) 



v„ k 



r k = In 



w b - w Vk w a - w Uk 
Wb - w Uk w a - w, 



(3.9) 



(3.10) 



(3.11) 



for k = l,...,g. Moreover, from the integral representation of q2(a,b) and K\(a,b) (see (2.12) 

1 



and (2.13)), using (3.2) and ( |3.3| ) one gets 

Q2(a,b) 

K x (a,b) 



k'(w a )k'(w b )(w a - w b ) 2 ' 

1 1 k"(w a ) 

k'(w a )(Wb- W a ) 2k' 2 (w a )' 



Putting z = in (2.5 ) and taking the limit e — > leads to 

q x {a,b) = -q 2 (a,b), 



(3.12) 
(3.13) 

(3.14) 



due to the fact that the theta function tends to one and that its partial derivatives tend to zero. 

(3.15) 



In the same way, taking the limit e — > in (2.6) one gets 

K 2 (a,b) = 



Kx(a,b) 



4 Degenerate algebro-geometric solutions of n-NLS 



One way to construct solutions of (1.2) is first to solve its complexified version, a system of 2n 
equations of 2n dependent variables , ?/>* | , 



Kk=l 
' n 



dib* d 2 ib* 

-i~ot + -qJ- + 2 (£ Mtj ^ = 0, 3 = 1, 



(4.1) 
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where ipj(x,t) and vpj(x,t) are complex valued functions of the real variables x and t. This 
system reduces to the n-NLS s equation (1.2) under the reality conditions 



j - 



J 



, n. 



(4.2) 



Algebro-geometric solutions of the system (4.1 ) were obtained in [22] by the use of the degenerated 



versions (2.5) and (2.6) of Fay's identity; these solutions are given by: 



Theorem 4.1. Let lZ g be a compact Riemann surface of genus g > and let f be a meromorphic 
function of degree n + 1 on !Z g . Let z a G C be a non critical value of f, and consider the fiber 
f~ 1 (z a ) = {ai, . . . , a n +i} over z a . Choose the local parameters k aj (p) = f(p) — z a , for any point 
p 6 TZ g lying in a neighbourhood of Qj . Let DdC 9 and Aj ^ be arbitrary constants. Then the 



following functions {^j}™ =1 and {V^l are solutions of the system (4-1) 

e(z 



ipj(x,t) 



A 



D + r, 



e(z-D) 

q 2 (a n+1 ,aj) 9(Z 



exp {i (—Ej x + Fj t)} , 
D - 



A; 



e(z-D) 



exp {i (Ej x — Fjt)} . 



(4.3) 



Here G denotes the theta function (2.2) with zero characteristic, and Z = i V an+1 x + i W an+1 1, 
where vectors V an+1 and W Qn+1 are defined in (2.3). Moreover, Xj = f^ j ^uj, where lo is the 
vector of normalized holomorphic differentials, ana the scalars Ej , Fj are given by 



Ej = Ki(a n+ i,aj), 



Fj = K 2 (a n+ i,aj) - 2} j qi(a n+1 ,a k ) 



(4.4) 



k=l 



The scalars q%,Ki for i = 1,2 are defined in (2.1)-(2.10). 



The proof of this theorem is based on the following identity: 



n+l 



E v ^ = < 



(4.5) 



fe=i 



which is satisfied by the vectors V Qfc associated to the fiber / (z a ) = {oi, . . . ,a n _|_i} over z a . 



We shall use this relation to construct solutions of (1.2) in terms of elementary functions. 



Remark 4.1. The relationship between solutions of the Kadomtesv-Petviashvili (KP1) equation 
(generalization of the KdV equation to two spatial variables, see, for instance, |?3) and solutions 
of the multi-component NLS equation was investigated in [22j . This relationship implies that all 



solutions of equation (1.2) constructed in this paper provide also solutions of the KP1 equation 
as explained in |22j . 



In the next section, solutions of (4.1) in terms of elementary functions are derived from 



solutions (4.3) by degenerating the associated Riemann surface TZ„ into a Riemann surface of 



genus zero. Imposing reality conditions (4.2), by an appropriate choice of the parameters one 



gets special solutions of ( 1.2 ) such as multi-solitons and breathers. To the best of our knowledge, 



such an approach to multi-solitonic solutions of n-NLS s has not been studied before. Moreover, 
breather and rational breather solutions to the multi-component case are derived here for the 
first time. 
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4.1 Determinantal solutions of the complexified n-NLS equation 

Solutions of the complexified scalar NLS equation in terms of elementary functions were obtained 
in |7j, when the genus of the associated hyperelliptic spectral curve tends to zero. For specific 
choices of parameters, they get dark and bright multi-solitons of the NLS equation, as well 
as quasi-periodic modulations of the plane wave solutions previously constructed in |21j. A 
direct generalization of this approach to the multi-component case is not obvious, due to the 
complexity of the associated spectral curve. To bypass this problem and to construct spectral 



data associated to algebro-geometric solutions (4.3) in the limit when the genus tends to zero, 



we use the uniformization map between the degenerate Riemann surface and the sphere. Details 
of such a degeneration were presented in Section 3. 

Let us discuss solutions of n-NLS in genus zero. Consider the following meromorphic function 
f(w) on the sphere: 

n+l 

W - W ai 
W - w bi 



(4.6) 



where w a . ^ Wb k for all j, k, w a . ^ w ak for j ^ k, and a G C. Without loss of generality, 
put a = 1. This function is of degree n + l on the sphere, hence it represents a genus zero 
(n+l)-sheeted branched covering of CP 1 . Recall that a meromorphic function / on the sphere is 
called real if its zeros as well as its poles are real or pairwise conjugate. 

If not stated otherwise, the local parameter in a neighbourhood of a regular point w a (i.e. 
f'(w a ) 7^ 0) is chosen to be k a (w) = f(w) — f(w a ) for any w lying in a neighbourhood of w a . 



Solutions of the complexified system (4.1 ) associated to the meromorphic function / (4.6) on the 
sphere are given by: 

Proposit ion 4.1. Let j, k G N satisfy 1 < j < n and 1 < k < g. Let f be a meromorphic 
function (4-6) of degree n + 1 on the sphere, with complex zeros {u^}™^ 1 and complex poles 
{^fejlfci ■ Let d G C 9 and Aj ^ be arbitrary constants. Moreover, assume that w Uk ,w Vk G C 
satisfy 

f(w Uk ) = f(w Vk ). (4.7) 



Then the following functions are solutions of the complexified system (4-D 

, det(T iil 



det(Tj j0 " 



exp{i(-£jx + Fjt)}, 



q 2 (an+i,aj) det(Tj_i) 
Aj det(Tj i0 ) 



exp{i (Ej x — Fjt)}. 



(4.8) 



For (3 = —1,0, 1, Tj^ denotes the g x g matrix with entries where z J k = — d^ + firj^k- 

Here = i V an+1) k % + i W an+U k t, where the scalars V an+lt k and W an+1 ,k are defined in (3.9) and 
(3.10), and tj^ is defined in (3.11) with w a := w an+1 and w b ■ 



given by 



«n+l 

Wa 



The scalars Ej and Fj are 



Ej = Ki(a n+1 ,aj), Fj = - (ifi(a n+ i,qj)) 2 + 2 q 2 (a w +i , Qfc), 



fc=i 



where q2(a n +i,aj) and K\(a n -\-\,a,j) are defined in (3.12) and (3.13). 
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Proof. Consider solutions (4.3) associated to a Riemann surface lZ g of genus g, and assume 
f(a,i) = for any 1 < i < n + 1. Pinch all ,4-cycles of the associated Riemann surface lZ g into 
double points, as explained in Section 3. After desingularization, the meromorphic function / 
of degree n + 1 on lZ g becomes a meromorphic function of degree n + 1 on the sphere, given in 



general form by (4.6). In the limit considered here, the theta function tends to the determinantal 



form (3.7). Quantities defined on the degenerated surface and independent of the variables x and 



t were constructed in Section 3.3 and are given in (3.9 )-(3.15 ). Condition (4.7) follows from the 



fact that double points appearing after degeneration of lZ g are desingularized into two distinct 
points w Uk and w Vk having the same projection under the meromorphic function /. Note that 



equation (4.5) holds in the limit, since by (2.3) and (4.6) one has 



n+l 



1=1 



f{Wu k ) f(Wv k ) 



which by (4.7) equals zero for k = 1, . . . ,g. 



(4.9) 
□ 



Remark 4.2. Functions (4.8) give a family of solutions to the complexified multi-component 



NLS equation (4.1 ) depending on 3n + g + 2 complex parameters: w ai , for 1 < i < n + 1, dj. 
for 1 < k < g, and Aj for 1 < j < n. 



Remark 4.3. The following transformations leave equation ( |4.1[ ) invariant 

ipj(x,t) — >4>j ((3x + 2(3\t,[3 2 t) exp{-i (\x + X 2 t)} , 
tp*(x,t) — > P 2 %b* (/3x + 2/3At,/3 2 t) exp {i (A x + A 2 1) } 



(4.10) 



where A = /i f3 1 for any /i£C and any /3 7^ 0. Such a transformation may be useful to simplify 
the expressions in the obtained solutions and thus to facilitate the numerical implementation. 



4.2 Multi-solitonic solutions of n-NLS 



Imposing reality conditions (4.2) on the degenerate solutions (4.8) of the complexified system, 
one gets particular solutions of (1.2) such as dark and bright multi-solitons. Dark and bright 



solitons differ by the fact that the modulus of the first tends to a non zero constant and the 
modulus of the second tends to zero when the spatial variable tends to infinity. Such solutions 
were obtained in [7] for the one component case by degenerating algebro-geometric solutions, 
and describe elastic collisions between solitons. Elastic means that the solitons asymptotically 
retain their shape and speed after interaction. The interaction of vector solitons is more complex 
than the one of scalar solitons because inelastic collisions can appear in all components of one 
solution (see for instance |T]). 



In what follows N £N with N >1. 



4.2.1 Dark multi-solitons of n-NLS s , s ^ (1, . . . , 1). 

Dark multi-soliton solutions of 2-NLS 5 were investigated in |35j . The dark iV-soliton solution 
derived here corresponds to elastic interactions between iV dark solitons. Moreover, it is shown 
that this type of solutions does not exist for the focusing multi-component nonlinear Schrodinger 
equation, i.e., in the case where s = (!,...,!). 
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Proposi tion 4.2. Let j,k 6 N satisfy 1 < j < n and 1 < k < N. Let f be a real meromorphic 
function (4-6) of degree n + 1 on the sphere, having n + 1 real zeros 

{wa^i Choose 9 G K and 



d G M. N . Moreover, assume that w Uk ,w Vk G C satisfy (4-1) and 



u ' 



W 



Vk ■ 



(4.11) 



Put Sj = sign(f'(wa n+l )f'(w a:j )). Then the following functions define smooth dark N-soliton 
solutions of n-NLS s , where s = (s±, . . . , s n ) with (!,...,!), 



ipj(x,t) 



(4.12) 



Here A 



/) I 1 / 2 , a^c? i/ie remaining notation is as in Proposition 4-1 with g = N . 



Proof. Let us check that the functions ijjj and ij)*- defined in (4.8) satisfy reality conditions (4.2) 
with Sj = sign(f'(w an+1 )f'(w a:j ))- Put Aj = |g 2 (a n +l, aj)\V*m (4.8). Then with the above 
assumptions, it is straightforward to see that ipj = Sj ipj where Sj = sign((72(a n +i, aj)), which by 



3.12) leads to Sj = s\gn(f'(w an+1 )f'(w aj )). Moreover, with (4.11) it can be seen that condition 



(4.5) is equivalent to 



W 



+ ^/'K +1 ) 



f'(w a .) \w a . 



w 



(4.13) 



for k = 1, . . . , N. Therefore, by (4.13) the quantity f'(w an+1 )f'(w a ) cannot be positive for all j, 



which yields s / (1, . . . , 1). The solutions are smooth since the denominator in (4.12) is a finite 
sum of real exponentials. □ 



Remark 4.4. The dark iV-soliton solutions (4.12) depend on N + 1 real parameters d^,9 and 



a real meromorphic function / (4.6) defined by 2n + 2 real parameters. The solitons are dark 
since the modulus of the t/jj tends to Aj when 16R tends to infinity. 



Example 4.1. With the notation of Proposition 4-1 and 4-%> functions ijjj (4-12) are given for 
N = 1 by 

1 -j- e Zi-di+r 3 

ibj(x,t) = Aj ^ — j— 



^(-EjX+Fjt) 



4.2.2 Bright multi-solitons of n-NLS s . 

Bright multi-solitons of the NLS equation presented in [7] were obtained by collapsing all branch 
cuts of the underlying hyperelliptic curve of the algebro-geometric solutions. This way they get 
solutions expressed as the quotient of a finite sum of exponentials similar to dark multi-solitons, 
except that the modulus of the solutions tends to zero instead of a non-zero constant when the 
spatial variable tends to infinity. Following this approach, a family of bright multi-solitons of 



n-NLS s is obtained here by further degeneration of (4.8). 



For the multi-component case there exist two sorts of bright soliton interactions: elastic 
or inelastic. Inelastic collisions between bright solitons were investigated in |36| for the two 
component case and in |23] for the multi-component case. The family of bright multi-solitons of 
n-NLS s obtained here describes the standard elastic collision with phase shift. Notice that there 
exist various ways to degenerate algebro-geometric solutions. Therefore, it appears possible that 
bright solitons with inelastic collision can be obtained by different degenerations. 
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Proposition 4.3. Let j G N satisfy 1 < j < n. Take w a ., 9 G R and choose d G C suc/i i/iai 
^2fe-i = ^2fc- Moreover, let w U2k ,w V2kl G C satisfy 



(4.14) 

/or 1 < /c < iV. Choose jj G K and put Sj = sign^j). Then the following functions give bright 
N-soliton solutions of n-NLS s 



ipj(x, t) = Aj e 



i6 



det(K i 



det(M) ' 

where Aj = It?] 1 / 2 |w aj | _1 . Here ~Kj and M are 2N x 2N matrices with entries (Kj)ik and 
given by: 

- for i and k even: (Kj) ik = 5 2 ,i ^ e^+^a+^k) + S _ § 2 . §2 k 



(4.15) 



Ilk 



w 



II k 



+ h,k ($2,: 



1) ^ e H z i- z ^+ f i,i-fj,2) 



Here a 



- for i even and k odd: 

- for i odd and k even: 

- for i and k odd: 

for i, k even, or i,k odd: 

- for i even and k odd: 

- for i odd and k even: 
where 



w„ 



w 
a, 



"2 



a 



II k 



K 



w 



■jjik — 

W Vi - W v 



^1 p \(zi+Zk+r 3 ,i+rj, k ) 



W 



a„. — a 



i'k 



w 



w Vi — w 



uk 



a 



U2k-1 



7 = 1 V 



w aj - w V2k _ x 



(4.16) 



Moreover, z k is a linear function of the variables x and t satisfying z 2 k = Z2k-\, given by 

z 2 k-i = ia^.i x + ial 2k _ i t - d 2k -\. 
The scalars fjf. satisfy rj2k = ~ r j,2k-i where 

fj,2k-i = In 



W n . A - W. 



V2k-1 



Waj W V2k _ 1 Oiu2k-i 



Proof. Consider functions (4.8) obtained from (4.3) for the choice of local parameters k ai : 



k ai {w) = { ll f'{w ai ))- 1 f{w) 
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for any w lying in a neighbourhood of w ai , i = 1, 



(4.5) becomes 



n+1 



n + 1, and assume g = 27V". Hence condition 
1 



w 







(4.17) 



for A; = 1, . . . , N, Now put Aj = \q2(a n +i, aj)! 1 / 2 in (4.8), where 

g 2 (a„+i, aj) = 7„+i 7j (w a „ +1 - ^a,) -2 - 
Choose a small parameter e > and define = — In e + cf^, for k = 1, . . . , 2iV, and 

w an+1 + e 2 a-i, 



«2fc-l 



^a n+1 +e 2 a U2 1 fc „ 1 , 



(4.18) 



for k = 1, ... , iV. Now put 7 n +i = e 2 and consider in the determinant det(Tj i i) appearing in 



(4.8) the substitution 



L 



2i 



L 



2i 



(^•,1)2,2 

for i = 2,...,N, where denotes the line number k of the matrix Tj i, and (Tj i)j k denotes 
the entries of this matrix. In the limit e — > 0, it can be seen that the functions ifjj given in (4.8) 



converge towards functions (4.15), where the following change of parameters (eliminating the 
parameter w an+1 ) has been made: 



Wa 3 + W an+1 , W U2k -»■ W U2k + W an+l , W V2k ^ -> W V2k _ x + W an+1 , 



for j = 1, . . . ,n and A; = 1, . . . ,N. Analogous statements can be made for the functions ip*. 
By assumption, it is straightforward to see that the functions ipj and ip* obtained in the limit 
considered here satisfy the reality conditions ipj = Sjipj with Sj = sign(7j). Moreover, in this 
limit condition (4.17) yields (4.16). □ 



Remark 4.5. The bright iV-soliton solutions (4.15) depend on 2N complex parameters efofc-i) 



w 



V2k- 



1 , and 2n + 1 real parameters w aj ,^fj,9. Moreover, all parameters appearing in (4.15) are 



free, contrary to the dark multi-solitons (4.12) where parameters w Uk and w Vk have to satisfy 
the polynomial equation (4.7). The solitons are bright since the modulus of the ipj tends to zero 



when i£R tends to infinity, in contrast to the dark solitons. 

4.3 Breather and rational breather solutions of n-NLS 

Solutions obtained here differ from the dark multi-solitons studied in Section 4.2.1 by the reality 



condition imposed on parameters w Ui and w Vi in solutions (4.8) of the complexified system for 
i = l,...,g. By an appropriate choice of parameters, one gets periodic solutions (breathers) 
as well as rational solutions (rational breathers). The name 'breather' reflects the behavior of 
the profile of the solution which is periodic in time (respectively, space) and localized in space 
(respectively, time) . This appears to be the first time that explicit breather and rational breather 
solutions of n-NLS s are given. 



In what follows N £ N with N > 1. 
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4.3.1 Multi-Breathers of n-NLS s . 



Multi-breather solutions of n-NLS s are given in the following proposition. The iV-breather solu- 
tion corresponds to an elastic interaction between N breathers. 

Proposi tion 4.4. Let j,k 6 N satisfy 1 < j < n and 1 < k < N. Let f be a real meromorphic 
function (4-6) of degree n + 1 on the sphere, having n + 1 real zeros {Wcn}?=i- Choose 6 6 WL and 
take d G C 2N such that d 



2fc-l 



d 2 k- Let Wu 2k i'Wu2k-ii w v2ki' w v2k-i G C satisfy (4.7) and 



w 



"2k 



w 



V2k-1 ' 



w 



V2k- 



(4.19) 



Put s 
n-NLS : 



j — sign(f'(wa n+1 )f'(w aj ))- Then the following functions define N -breather solutions of 



^OT exp{i( -^ + Fjt>} ' 



(4.20) 



where Aj = \q2(a n +i, a-,-)! 1 / 2 , and the remaining notation is the same as in Proposition 4-1 f or 



9 



2N. 



Remark 4.6. Functions (4.20) cover a family of breather solutions of n-NLS s depending on N 
complex parameters dk, a real parameter 6, and a real meromorphic function / (4.6) defined by 
2n + 2 real parameters. 



To simplify the computation of the solutions, we apply transformation (4.10) to the solutions 

1 



(4.20), with (3 and A given by 

P 



1, 



A = 7: f" (W an+1 ) f {W an+1 )' 



(4.21) 



Hence, the quantity /"(w a „ +1 )/'(w«„+i 

2- 



Va n+1 ,k i n the expression (3.10) for the scalar W an+lt k 



disappears, as well as the quantity hf"(w an+1 )f'(w an+1 ) 2 in the expression (3.13) for the scalar 



Ki(a n+ i,a,j). 

Example 4.2. Figure^shows a breather solution of the 4-NLS s equation with s = (— , — ,+, — ). 
It corresponds to the following choice of parameters: w ai = 10, w a2 = —5, w a3 = —1/3, w ai = 
1/4, w a5 = 1/2, andw ui w 0.55 — 0. Hi with f(w Ul ) = 2i, w U2 » — 0.35+0. 07i with f(w U2 ) = — 2i. 

Example 4.3. Figure shows an elastic collision between two breather solutions of the 4~NLS S 
equation with s = (—,+,+,—). It corresponds to the following choice of parameters: w ai = 



1/3, w 



"2 



3, w 



1/7, W. 



a 4 



2, W a5 = !> w bi 



-1, w b2 = 4, w b3 



-2, w b . 



w„ 



0.55 — 0. Hi with f(w Ul ) = 2i 7 w h 



-0.35 + 0.07i with f(w U2 ) = —2i, w v 



0, and 
-0.91-0.52i 



with f{w. 



U3 ) 



10 — 5i, and w Ui ~ 14.46 + 5.32i with /(m 



U4 i 



10 + 5L 



4.3.2 iV-rational breathers of n-NLS s , for 1 < N < n. 



Here we are interested in solutions of n-NLS s that can be expressed in the form of a ratio of 
two polynomials (modulo an exponential factor). These solutions, called rational breathers, are 
neither periodic in time nor in space, but are isolated in time and space. They are obtained 
from breather solutions (4.20) in the limit when the parameters w V2k _ l and w u ^,, , tend to each 



J u 2k -i 



others, as well as the parameters w V2k and w U2k , 



for k = 1, . . . N. An appropriate choice of the 
, 2N, leads to limits of the form 0/0 in the expression for the 
breather solutions. Thus, by performing a Taylor expansion of the numerator and denominator 



parameters di in (4.20) for i = 1 



in (4.20), one gets a family of iV-rational breather solutions of n-NLS s . 
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Proposi tion 4.5. Let N,j E N satisfy 1 < N < n and 1 < j < n. Let f be a real meromorphic 
function (4-6) of degree n + 1 on the sphere, having n + 1 real zeros {waj^i Choose 9 eRand 



take d £ C 2N such that d 2 k = d 2k -i f or 1 < k < N . Moreover, let w U2k _ 1 ,w V2k EC, 1 < k < N, 
be complex conjugate critical points of the meromorphic function f , i.e., they are solutions of 
f'{w) = 0, which is equivalent to 



n+l 



E 7 ^2 = °- ( 4 - 22 ) 

Put Sj = sign(f'(w an+1 )f'(w a )). Then the following functions give N -rational breathers of n- 
NLS S 

Mx, t) = Aj e ie exp {i (-Ej x + Fj t)} , (4.23) 



where Aj = \q 2 (a n +i, aj)\ 1 / 2 . For /3 = 0, 1, Kj j( g denotes a 2N x 2N matrix with entries (Kj p)i jk 
given by: 

- for i and k even: {Kj^) ik = (1 - 5 ijk ) k,k {z k + /3fj, k ) 



w Vi - w Vk 



1 



- for i even and k odd: (Kj^)i k 

- for i odd and k even: (Kj,p)ik 

- for i and k odd: (K,- p) ik = - (1 - 6 i k ) 6 i k (z k + f3 f j k ). 

w Ui ~ w Uk 

Here z k is a linear function of the variables x and t given by 

Zk = i Va n+1 ,k x + i W an+1 ,k t - d k 



for k = 1, ... , 2N, where V an+lj2k = ~ K n+1 ,2fc-i and W an+lj2 k = - W an+lj2k -i with 
1 1 12 



for k = 1, . . . , N. Scalars fj^ k satisfy Tj t2k = — 2fc— l an d are given by 



W an+1 ~ W aj 



(w an+1 ~ Wu^) (w aj - Wu^) 
Scalars Ej,Fj are defined by 

1 n 
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Proof. To symplify the expression for the obtained solutions, apply the transformation (4.10) 
to functions (4.20) with (3 and A as in (4.21). Let e > be a small parameter and define 

+ <"„ J; . (4.24) 



dk = edk + 17T, for k = 1, . . . , 2N. Moreover, assume 



w 



V2k-1 



w 



"2fe- 



.i + ea 



V2k-i: 



W 



U2k 



W 



'-2k 



for some ct V2k _ i: a U2k £ C, where k = 1, . . . , N. Note that equation number k of system (4.9) 
can be written as 



n+1 



f(w Vk )f(w Uk ) 



pi f'(Wa 3 ) (Wa 3 - W Vk ) (W aj - W Uk ) 



f{Wv k ) ~ f(Wu k ) 



W 



w 



(4.25) 



y-k 



Hence, in the limit e — > 0, equation (4.25) becomes 

n+1 



y— 
^ fid, 



V2k-1/ 



i f'{w aj ) (w aj -w V2k _J< 



V2k-ih 



and 



n+1 



y— 

pif'(w aj ) K 



10, 



fiw. 



U2k)' 



V-2kJ 



for k = 1, . . . , iV. Therefore, choose w V2k _ 1 and ti^k to be distinct critical points of the meromor- 
phic function / for k = 1, . . . N, i.e., they are solutions of f'(w) = 0, in such way that equation 
(4.5) holds in the limit considered here. Since the condition f'(w) = is equivalent to solve a 
polynomial equation of degree 2n, it follows that 1 < N < n. Now take the limit e — > in (4.20). 
Note that parameters o. V2k _ l , o. U2k cancel in this limit, and the degenerated functions take the 
form ||03i. □ 



Remark 4.7. Functions (4.23) provide a family of rational breather solutions of n-NLS s de- 
pending on N complex parameters dk, a real parameter 6, and a real meromorphic function / 
(4.6) defined by 2n + 2 real parameters, chosen such that / admits complex conjugate critical 
points. 



Example 4.4. With the notation of Proposition 4-5 the functions ipj (4-23) for N = 1 are given 
by 



exp {i (—Ejx + Njt)} , 



where B = (2Im(w Ul )) 



-2 



Example 4.5. Figure [3] shows a rational breather solution of the 4-NLS 8 equation with s = 
(+,+,+,+). It corresponds to the following choice of parameters: k ak (w) = f'(w a .)f(w) for 
k = 1, . . . , n + 1, with w ai = 3, w a2 = 5, w a3 = 7, w a4 = 0, w a5 = 4, and w Ul « 4.53 + 0.56i 
being a solution of Y^i=ii w ~ w cn)~ 2 = 0- We observe that functions ip2 an d ip3 coincide with 
the Peregrine breather well known in the scalar case [33], whereas functions belong to a 

new class of rational breathers which does not exist in the scalar case. This new type of rational 
breathers emerges due to the higher degree of the meromorphic function associated to the solutions 
of n-NLS s for n> 1 . 



18 




19 



Example 4.6. Figure shows a 2-rational breather solution of the 4~NLS S equation with s = 
(+,+,+,+). It corresponds to the following choice of parameters: k ak (w) = f'(w ak )f(w) for 
k = 1, . . . ,n + 1, with w ai = 3, w a2 = 5, w a3 = 7, w ai = 0, w a5 = 4, and w Ul « 4.53 + 0.56i, 
w uz ~ 3.45 + 0.56i being solutions of Y27=i( w ~ w ai)~ 2 = 0, and dk = 10. Variation of the 
parameters dk leads to a displacement in the (x,t)-plane of the rational breathers appearing in 
each of the pictures of Figure [^} 



5 Degenerate algebro-geometric solutions of the DS equations 



Solutions of the DS equations ( 1.3 ) in terms of elementary functions constructed here are obtained 



analogously to the solutions of the n-NLS equation, therefore some details will be omitted. Let 
us introduce the function <fi := <I> + /o| -i/^ | 2 , where p = ±1, and the characteristic coordinates 
^ = \{x — iay), rj = ^(x + iay), a = i, 1 in (1.3). In these coordinates the DS equations become 



1 



0. 



(5.1) 



Recall that DSP denotes the case a = i (here £ and ry are both real), and DS2 P the case a = 1 
(here £ and r\ are pairwise conjugate). 



To construct solutions of (5.1) in terms of elementary functions, let us first introduce its 
complexified version: 



i i>t + 2 + Vv?) + 2 <P V> = 0) 
iV<t* + ^ + ^) + 2^*=0, 



(5.2) 



where 99 := $ + ^^*. This system reduces to (5.1) under the reality condition: 



ip* = pip, 



(5.3) 



which leads to tp = (j). Theta-functional solutions of (5.2) were studied in [22] and can be written 
in the following form. 

Theorem 5.1. Let 7Z g be a compact Riemann surface of genus g > 0, and let a,b £ 7Z g be 
distinct points. Take arbitrary constants DgP and A, K\,K2 € C \ {0}, h € C. Denote by £ a 
contour connecting a and b which does not intersect cycles of the canonical homology basis. Then 



for any £, r/, t S C, the following functions ip, if)* and <p are solutions of system 1 5.2) 

6(Z - D + r 



Vi^t): 



0(Z-D 

k\k 2 92(0, b) 6(Z - D - r 



A 



exp{-i(G 1 e + G 2 r ? -G 3 §)} ) 

exp{i(Gi£ + G 2 7/-G 3 |)} 



e(z-D) 



(5.4) 



^,77,*) = ^ (ln6(Z - D))« + - (ln9(Z - D)) w + J. 
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Here Z = iKiV a £ — i^V^ + W a — tt^Wb) 55 where the vectors V a , V;, and W a , 



were introduced in (2.3). Moreover r = J f u, where u is the vector of normalized holomorphic 



differentials, and the scalars G\,G 2 ,G 3 are given by 



G\ = ki Ki(a, b), G 2 = n 2 Ki(b,a) 
G 3 = k\ K 2 {a, b) + k\ K 2 (b, a) + h. 



(5.5) 
(5.6) 



Scalars q 2 (a,b), Ki(a,b), K 2 (a,b) are defined in (2.8), {2.ty , (2.10) respectively. 

Remark 5.1. In the case where vectors V a and satisfy V a + = 0, as mentioned in |22j . 



solutions of the Davey-Stewartson equation become solutions of the NLS equation (1.1) under 
an appropriate change of variables. 



In this section, we study the behaviour of theta-functional solutions (5.4) of the complexified 



DS equations when the Riemann surface degenerates into a Riemann surface of genus zero. 



Imposing the reality condition (5.3), for particular choices of the parameters one gets well-known 



solutions such as multi-soliton, breather, rational breather, dromion and lump. This appears to 
be the first time that such solutions of DS are derived from algebro-geometric solutions. 

5.1 Determinantal solutions of the complexified DS equations 



Here solutions of the complexified system (5.2) are given as a quotient of two determinants. In 



the next subsections, this particular form will be more convenient to produce special solutions 



of the DS equations (5.1). 



Proposition 5.1. Let k £ N satisfy 1 < k < g. Let w a ,w b ,w Uk ,w Vk ,h £ C, and A,ki,k 2 £ 



C \ {0}. Choose d G C g . Then the following functions are solutions of the system (5.2) 



^,V,t)=A^^exp{-i(G 1 Z + G 2V -G 3 l)} 



Kl K 2 



det(T_ 



A{w a -w b ) 2 det(To) 



^ exp{i(G 1 t; + G 2 r ) -G 3 l)} 



(5.7) 



1 1 h 

<p(£,ri,t) = o (lndet(T )) K + r (lndet(T )) w + -. 



the scalars r^ are given in (3.11) and 



For (5 = —1, 0, 1, Tg denotes the g x g matrix with entries (3.8) with zj~ = — dk + Pr^. Here 

(5.8) 



Z = i Kl V a £ - i K 2 V b T] + i (k? W a - k\ W b 



with 



1 



1 



c,k 



W r — VJ. 



W r — W. 



w, 



c,k 



+ 



1 



(w c -w Vk ) 2 (w c -w Uk ) 



where c G {a, b}. The scalars G\,G 2 , G 3 are given by 

Kl „ k 2 



Go 



G, 



w b - w a 



W a ~ W b 



G\-G 2 2 + h. 



(5.9) 



(5.10) 
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Proof. Consider solutions (5.4) of system (5.2) in the limit when the Riemann surface degener- 



ates to a Riemann surface of genus zero, as explained in Section 3. In this limit, choose the local 
parameters k a and kf, near a £ TZq and b G TZq to be the uniformization map between the degen- 
erate Riemann surface IZo and the w-sphere. Hence, for any w £ TZq lying in a neighbourhood of 
w a £ IZq, k a (w) = w — w a - Therefore, quantities independent of variables £, r/ and t are obtained 
from OM3l5i. 



□ 



Remark 5.2. Functions (5.7) give a family of solutions of the complexified system, involving ele- 
mentary functions only. These solutions depend on 3g + 6 complex parameters w a , Wb, h, A, k±, K2 
and w Uk , w Vk , dk- Varying these parameters we will obtain different types of physically interesting 
solutions investigated in the next subsections. 



5.2 Multi-solitonic solutions of the DS equations 

Soliton solutions of the DS equations were shown to be representable in terms of Wronskian 
determinants in |5]. Single soliton and multi-soliton solutions corresponding to the known one- 
dimensional solutions can be obtained from this representation. These solitons are pseudo-one- 
dimensional in the sense that in the (x, y)-plane, they have the same form as one-dimensional 
solitons in the (x, £)-plane, but that they move with an angle with respect to the axes. The 
multi-soliton solution describes the interaction of many such solitons each propagating in differ- 
ent directions. 



In what follows N £ N with N > 1. 



5.2.1 Dark multi-soliton of BS1 P and DS2 



Here dark multi-solitons of the DS1 P and DS2 + equations are derived from functions (5.7) for 
an appropriate choice of the parameters. They were investigated in |41j . 



Put g = N and A = \k,iK2\ 1 ^ 2 \w a — io&| 1 in (5.7). Moreover, assume ftgi and d £ 



Reality condition for DS1 P . Let us check that with the following choice of parameters, 

w a ,w b £R, Ki,K 2 £l\{0}, w^ = w Uk , k = l,...,N, (5.11) 



functions ip and ijj* in (5.7) satisfy the reality condition ip* = pxfj with p = — sign(«;iAC2)- Indeed, 
this can be deduced from the fact that G\, G2, G3 £ M, and 



det (Tjg) = det (l> ) = det (T_ 



(5.12) 



since u and v can be interchanged in the proof of ( |3.7[ ). Therefore, functions ip and <f> in ( |5.7| 
define dark multi-soliton solutions of DS1 P . 



Smoothness. The dark multi-soliton solutions obtained here are smooth because the denomi- 



nator det (To) of functions if) and <f> (5.7) consists of a finite sum of real exponentials (see ( |3.7[ )), 
since £, 77, t are real. 

Remark 5.3. One gets a family of smooth dark multi-soliton of the DS1 P equation, depending 
on N + 6 real parameters w a , io&, h, K%, K2, dk, a phase 8, and N complex parameters w Uk . 
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Reality condition for DS2 + . Let us check that with the following choice of parameters, 



Wb, 



w u ,,w v , G 



k = 1, 



(5.13) 



the functions ?/> and (5.7) satisfy the reality condition ip* = ip. With (5.13), it is straightfor- 



ward to see that (5.12) is als o satisfied. Moreover, since d = G 2l G 3 G R and (w a - w h y < 0, 



the functions ip and ip* (5.7) satisfy the reality condition ip* = ip. Therefore, they define dark 



multi-soliton solutions of DS2 . 

Smoothness. To get smooth solutions, additional conditions are needed to ensure that det(To) 
does not vanish for all complex conjugate £ = fj. For instance, if 



w Vl < w Ul < w V2 < w U2 < 



< w VN < w UN , 



the scalars (B)^ (3.6) are real for any i, k G {1, . . . , N}. Therefore, the functions ip and (p (5.7) 
are smooth, since their denominator does not vanish as a finite sum of real exponentials. 

Remark 5.4. One gets a family of smooth dark multi-soliton of the DS2 + equation, depending 
on 3iV + 1 real parameters h, w Uk , w Vk , d^, a phase 9, and 2 complex parameters w a ,K\. 

5.2.2 Bright multi-soliton of DSl'' and DS2~ 

In this part we construct bright multi-soliton to the DS1 P and DS2 - equations. It is well known 
that such solutions can be written in terms of a quotient of sums of exponentials, for which the 
modulus tends to zero if the spatial variables tend to infinity. 



To get bright multi-soliton solutions, one degenerates once m ore solutions (5.7) of the com- 
plexified system. Put g = 2N and A = \k\K2\ 1 ^ 2 \ w a — Wb\~ l in (5.7), and take li£M. 

-me + 4, for k = 1, . . . , 2N, 



Degeneration. Choose a small parameter e > and define d^ 
and 



w 



Wa +6%^ (w a ~ Wb), 



W 



V2k-1 



Wb + ea 



{w a - w b ), 



w 



II '2k 



[w„ 



w b ), 



w 



I'M 



w a + e a* {w 



(5.14) 



w b ) 



for k = 1, . . . , N. Moreover, put K\ = ek\ (w a — Wb), and k 2 = ek 2 (w a — Wf,). Consider in the 

L 2 



determinant det(Ti) appearing in (5.7) the substitution 

(Ti) 2i ,2 



L 



2 1 



L 



2i 



(Ti) 



2.2 



for i = 2, . . . , N, where denotes the line number k of the matrix Ti and (Ti)j k the entries of 
this matrix. An analogous transformation has to be considered for the matrix T_i appearing in 



function ip*. Now take the limit e — > in (5.7). The function ip obtained in this limit has the 



form (5.17). Notice that in this limit, the dependence on the parameters w a and Wb disappears. 



Reality condition for DS1 P . It is straightforward to see that, with the following choice 
of parameters, 



2k- 



a 



«2fe-l 



a 



V2k ' 



a. 



a 



V2k-1 ' 



A: 



1, 



N, (5.15) 
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the functions ip and tp* obtained in the limit considered here satisfy the reality condition ip* = pip 
with p = — sign(£iK2)- 



Reality condition for DS2 . In the same way, with the following choice of parameters 
~~ k = 1 



k\ — K2-, d,2k-l — ^2fc* a u 2k -! — a 



«2fc' 



a 



V2k-1 



a 



5***5 5 



the functions ip and ip* obtained in the considered limit satisfy the reality condition ip* 



(5.16) 



The solutions. Let 9 E M. With (5.15), the following functions of the variables £,77, £ ob- 



tained in the considered limit, give bright iV-soliton solutions of DS1 P where p = — signal £2) 
and 7 = 0; because of (5.16) these functions define bright ./V-soliton solutions of DS2~ where 
7 = 1: 



w det(K) 
det(M) ' 



0(e,r/,t) = i(lndet(M))^ + ^(lndet(M)) w + ^ 



(5.17) 



where A = |kiK2| 1//2 - Here K and M are 2N x 2iV matrices with entries (IC)jfc and (M)jfc given 
by: 



- for i and k even: (K) 



+ $2,k (S24 - l)e 

„2 



5 (zi- 22+^-7*2) 



- for i even and k odd: 

- for i odd and A; even: 

- for i and A; odd: 

- for i, k even, or i, k odd: 
- otherwise: 



K) ifc 
K) ifc 



n 



u 7 . 



a. 



°^ e |(^+2 fc +n+r fe ) 



at 



Hit 



-1) 



+1 a ^ a u k e \{zi+z k ) 



a, 



a 



Here ^ is a linear function of the variables £, rj and i given by 

S.2 „,2 



2 



17T 



^2fc = - i % at> 2fc £ - i«2 a« 2fc ?7 - i («i + &2 «u 2 J 7, - &2k + 7 
for A; = 1, . . . , N . Moreover, the scalars are defined by 

f k = {-l) k ln{-a Vk a Uk } , & = !,..., 2N. 



ITT 



Remark 5.5. i) With (5.15 ), functions (5.17) give a family of bright multi-soliton solutions of the 



DS1 P equation depending on 3iV complex parameters d,2k-i,o: U2k _ 1 ,a U2k and 4 real parameters 

h,9,Kl,K2- 



ii) With (5.16), functions (5.17) provide a family of bright multi-soliton solutions of the DS2 



equation depending on 3N+1 complex parameters c^fc-i, a u 2 k-i > a v 2 k-i • -^l anc - 2 rea l parameters 
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5.3 Breather and rational breather solutions of the DS equations 



The breather solutions of the DS equation were found in |40| . Here a family of breather solutions 
and rational breather solutions of the DS1 equation are derived from algebro-geometric solutions. 
These solutions resemble their 1 + 1 dimensional analogues. In particular, the profiles of the 
corresponding solutions of the DS equation in the (x, y, t) coordinates look as those in the (x, t) 
coordinates extended along a spatial variable y. 



5.3.1 Multi-Breathers of DSl^ 



The iV-breather solution obtained here corresponds to an elastic interaction between ./V breathers. 



Put g = 2N and A 



\KlK2 



1/2 



u>b\ 1 in (5.7). It is straightforward to see that with the 



following choice of parameters, 

W a , Wb, /l£l, Kl,K2 Sl\{0} 



l 2k-l 



l 2k, 



W 



V2k 



w 



"2fe-l ' 



w 



V2k-1 



w 



«2fc' 



for k = 1,...,N, functions tp and tp* (5.7) satisfy the reality condition ip* 



(5.18) 
pip with p = 



— sign(«iK2)- Therefore, analogously to the n-NLS equation, functions ij) and <f> in (5.7) give 
iV-breather solutions of DS1 P . 

Remark 5.6. One gets a family of breather solutions of DS1 P depending on 3iV complex pa- 



rameters d,2k-l, w % 



and 6 real parameters w a , Wb, h, K\, H2 and a phase 9. 



Example 5.1. Figure^ shows the evolution in time of 2-breather solution of DS1 with the 



following choice of parameters: w a 
1 + 4i, Ki = K2 = 1, dk = h = 0. 



i, w b 



-1, w. 



"1 



5 - 2i, w U2 



2 + i, w. 



1, W Ui 




Figure 5: 2-breather of DS1~ at a) t = 0, b) t = 45. 
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5.3.2 Multi-rational breathers of DSl^ 



In this part, we deal with rational solutions (modulo an exponential factor) of the DS1 P equation. 
These solutions are obtained as limiting cases of the breather solutions. The iV-rational solutions 
describe elastic interaction between N rational breathers, and are expressed as a quotient of two 
polynomials of degree N in the variables £, 77, t. 



Assume g = 2N and put A = {kik^ 1 ^ 2 \w a — 1 in (5.7). 

Degeneration. Let e > be a small parameter and define d k = ed k + ivr, for k = 1, . . . , 2N, 
and 

w v 2k -i = w « 2 ft-i + e a v 2 k-i 5 w U2k = w V2k + e a U2k (5.19) 
for k = 1, . . . , N . It is straightforward to see that det(T^) « e 2Ar Pg, where Pp is a polynomial 
of degree 2N with respect to the variables £,77 and t. Now take the limit e — )• in (5.7). The 
function if) obtained in this limit is an iV-rational breather solution of DS1 P given by (5.21). 

Reality condition. Imposing the following constraints on the parameters: 

w a ,w b ,h£R, «i,k 2 eK\{0}, d 2k = d 2 k-l, w U2k _ 1 = w V2k , k = l,...,N, (5.20) 



it can be seen that the functions ip and \p* (5.7) in the considered limit satisfy the reality con- 
dition tp* = pip, with p = — sign(KiK2)- 

The solutions. Let 6 M. Then the following degenerated functions define iV-rational breather 
solutions of DS1 P 

V,t) = A e w exp { -i (G 1 £ + G 2 r, - G 3 |)} , 

H^r],t) = i(lndet(K ))^ + i(lndet(K )) w + ^ (5.21) 
where Kg, with j3 = 0, 1, is a 2iV x 2iV matrix with entries (Ka)^ given by 

- for i and fe even: (Kg) ife = (1 - 8 i>k ) S i)k (z k + f3 r k ) 



w Vi - w Vk 



1 



for i even and k odd: (ICg)^ 
for i odd and k even: (IC^)^ 

- for i and odd : (K^) ik = - (1 - 5 ijfc ) 5 ijk (z k + f3f k ). 



W U ; — W 



Here z k is a linear function of the variables £, 7/ and i given by 



Zfc = i KiV atk £ - i K 2 V b ,k V + i ( «1 Wa,Jfe ~ «2 W&,fc ) ^ ~ 4- 



Moreover, for c G {a, 6}, the scalars V c ^ k ,W C)k and r* satisfy V c>2k = V Cj2k -i, W C)2k = W c ^ k -i 
and f 2k = r 2k —i, and are given by: 

1 - 2 

Vc,2ifc-1 = 7 To") W Ci2k -i 
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r-zk-i 



W a ~ W b 



W 



"2fc- 



for fe = 1, ...,N. Constants G^G^Gs are given in (5.10). 



Remark 5.7. Functions (5.21) give a family of rational solutions of DS1 P depending on 2N 



complex parameters d,2k-i,w U2k _ 1 and 6 real parameters w a , h, 9, k%, k.2- 

Example 5.2. Figure shows the evolution in time of the 2-rational breather solution of DSF 
with the following choice of parameters: w a = 2, w\, = 1, w Ul 
1, dfe = h = 0. 



2i, io U3 = 2 + i, ki = k 2 



K 'HK -h+ 



Figure 6: 2-rational breather of DS1 at a) i = —5, b) t = 0, c) £ = 5. 

Example 5.3. Figure [?| (resp. Figure^! shows the interaction between a line rational breather 
and a rational breather solution of DS1~ with the following choice of parameters: w a = 2, Wb = 
—2, w Ul = 3i (resp. w Ul =3i + l) , w U3 = 2i, k\ = k<i = 1, df, = h = 0. By line rational breather 
we denote a growing and decaying mode localized only in one direction. 



5.4 Dromion and lump solutions of the DS equations 

Here we construct the dromion solution of DS1 P and the lump solution of DS2 - which correspond 
to solutions l ocal ized in all directions of the plane. These solutions arise by suitable degenerations 
of solutions (5.7) to the complexified system, and by imposing the reality condition ip* = pip. 
This appears to be the first time that such solutions are obtained as limiting cases of theta- 
functional solutions. 



5.4.1 Dromion of DSl^ 

Boiti et al. [8] have shown that the DS1 equation has solutions that decay exponentially in 
all directions. The solutions they obtained can move along any direction in the plane, and the 
only effect of their interactions is a shift in their position, independently of their relative initial 
position in the plane. Later, Fokas and Santini |15H39| pointed out that by an appropriate choice 
of the boundary conditions, the localized solitons (called "dromions") of the DS1 equation pos- 
sess properties which are different from the properties of one-dimensional solitons, namely, the 
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Figure 7: Interaction between a line rational breather and a rational breather of DS1~ at a) 
t = —50, b) t = —20, c) t = —5, d) t = 0, e) t = 10, f) t = 50. The rational breather propagates 
in the same direction as the line breather. 




Figure 8: Interaction between a line rational breather and a rational breather of DS1~ at a) 
t = —50, b) t = —20, c) t = —5, d) t = 0, e) t = 10, f) t = 50. The rational breather propagates 
transversally to the direction of the line breather. 
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performed solutions do not preserve their form upon interaction. For a particular choice of their 
spectral parameters, they recovered solutions previously derived by Boiti et al. For details on the 
theory of dromion solutions the reader is referred to |37j and references therein. In this section 
we explore how the simplest dromion solution can be derived from algebro-geometric solutions. 



Let us consider solutions of the complexified system obtained in (5.7). Assume g = 4 and 
put A = |kiK2| 1//2 \w a — Wbl' 1 - 



Degeneration. Choose a small parameter e > and define 
and 



ln(e) + df. for k = 1, . . . , 4, 



w 



ui 



ea. 



w 



w a + ea 



U2 > 



w 



«3 



Wb + ea v 



w 



U4 



ea,. 



w„ 



w a + e a vi , 



w r 



ea, 



e a 



V31 



w b + ea Vi . (5.22) 



Moreover, put K\ = ek\a Vl and K2 = ek2a U3 . Now consider the limit e — > in (5.7). The 



functions ip and <p obtained in this limit are given by (5.24). 



Reality condition. Choose w a ,w b , h,0 £l and k\, k% E K \ {0}. Moreover, assume 



dik — <hk-U a v 2k -i — a 



«2fc ' a t>2k 



a 



1,2. 



(5.23) 



Put p = — sign(£iK2). With (5.23), it can be seen that the degenerated functions ip and ip* 
obtained in the considered limit satisfy the reality condition ip* = pip. Therefore, the following 
degenerated functions give the dromion solution of DS1 P 



iP(£,V,t) = Ae 



<t>(£,V,t) = ldtfln{cp(£,r),t)} + J d m In V, *)} + T, 



(5.24) 



where 

<p(£,ri,t) = 1 + Ax e 2Re( Z1 ) +j 4 2e 2Re( Z3 ) + A3e 2Re( ^1 )+2Re( 2 3)_ 
Here Zk is a linear function of the variables £, rj, t given by 



K? t 



a Vl a?„ 2 



%3 



K2 k\ t 

cc,„ 2 



O: 



"3 



"3 



Constants A, A\, A2 and A3 are given by 



j4 = \k1k2 



1/2 



U7 a W 6 



(a 



V3 



Ao 



W b 



4Im(a„ 3 ) lm(a U3 ) 



au\ ) a Vl a U3 



A 3 = Ai A 2 + 



At 



w„ 



4Im(a^)Im(a Ml ) 
w a w b 1 



4Im(a„ 1 ) lm(a U3 ) \ a%i\ a^ 3 1 



Moreover, in the case where A\ > 0, A2 > and A3 > 0, functions (5.24) are smooth solutions 
of DSP. 
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Remark 5.8. i) Functions (5.24) define a family of dromion solutions of DS1 P depending on 6 
complex parameters di, ds, a Ul , a Vl , a Ua , a V3 and 6 real parameters w a , wi, k\, k>2, h, 9. 
ii) In the case where a Ul ,a V3 G R, one gets localized breathers, namely, the solution oscillates 
with respect to the time variable (modulus of ip is constant with respect to t). 

Different degenerations can be investigated for larger values of g. The performed functions 
lead to particular solutions such as dromions which move along sets of straight and curved 
trajectories, as well as oscillating dromion solutions. We do not discuss these solutions here. 

5.4.2 Lump of DS2~ 

The lump solutions were discovered in |26j for the KP1 equation, and have been extensively 
studied. Arkadiev et al. [6j have constructed a family of travelling waves (the lump solutions) 
of DS2~ that we rediscover here. 



Let us consider functions \fj,ip*,4> given in (5.7), assume g = 2 and put A = \k\K2\ 1 ^ 2 



\w n 



Wb\ 1 - Moreover, consider the following transformation which leaves the system (5.2) invariant: 



^(£, T i ,t)^nl>(Z + p 1 t,ri + ft t, t) exp { -i (ft i + ft V + [Pi + /3 2 2 ) |) } , 
r^t)^^* (£ + ft^ + ft M ) exp {i(ft£ + ft77+(ft 2 + ft 2 ) §)}, 

V, t) -»• <j> (£ + ft t, T) + ft t, t) , (5.25) 

where ft = U{ k~ 1 for some /Xj G C. 

Degeneration. Choose a small parameter e > and define dk = ivr + edf,, for k = 1,2, 
and 

Wvi = w a + e a Vl , w Ul = w a + ea Ul , 

w V2 = w b + ea V2 , w U2 = w b + e a U2 . 



Moreover, put = e k^, and //j. = e jik for = 1,2. Now take the limit e — >• in (5.25). The 



functions ip and obtained in this limit are given in (5.26). 



Reality condition. Choose w a ,Wb G C such that w a = — Wb or w a ,Wb G M. Take h, G R and 

assume 

«1 = «2, fii = £12, d\ = d2, = a V2 , a Ul = a U2 . 

With this choice of parameters, it can be seen that the functions ip and ip* obtained in the limit 
considered here satisfy the reality condition i/j* = — ip. 

The solutions. Therefore, the following degenerated functions provide smooth solutions of 
DS2~ 

Ae ie 

4(x, y, t) = exp { -i (2 Re(ft + Re(ft 2 ) t) } , 

B + \zi\ A 

0(x,y,t) = ^d^]n{B + \z 1 \ 2 } + ^d^\n{B + \z 1 \ 2 } + ^, (5.26) 
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where £ = x + iy and B\ = \i\ /% . Here z\ = i V a ,\ {k\ i + fi-it) — d\ and 

1*11 1 2 I 

va,i = > ^ = i 1 ■ ts - 



a ui a vi \w a - Wb\ a Vl a Ul (wb - w a )* 



Simplifications. To simplify (5.26), put 



j 1 Ol ul Q! vl 

d 1 = - 7 r — , V = — — r, \ = Pl, 

ValKl \Kl\\w a -W b \ 



for arbitrary /j, G C. In this way, functions (5.26) become 



exp{-2iRe(A£) -iRe(A 2 )t + W} 
ip(x, y,t) = u 12 i i 12 ' 

y, t) = ^ % In {|e + At + H 2 + W\ 2 } + \ %ln {|£ + A t + /x| 2 + M 2 } + (5.27) 



where £ = x + iy. Here A, ^, are arbitrary complex constants, and 6, h G M. Solutions (5.27) 
coincide with the lump solution previously obtained in [6]. 

6 Outlook 

In this paper, various classes of solutions to the multi-component NLS equation and the DS equa- 
tions in terms of elementary functions have been presented as limiting cases of algebro-geometric 
solutions discussed in a previous paper |22j . We did not construct all families of solutions present 
in the literature, but we believe that different degenerations will lead to interesting new or known 
solutions that are not presented here. 

In particular, future investigations might address bright multi-solitons of n-NLS with inelastic 
collision. This novel type of inelastic collision, which is not observed in 1 + 1 dimensional soliton 
systems, follows from a family of bright soliton solutions having more parameters than the ones 
presented here with standard elastic collision. We believe that also this kind of solutions arises 
from algebro-geometric solutions after suitable degenerations. 
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